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We attempt to reconcile large trilinear R-parity violating interactions in a supersymmetric (SUSY) 
theory with the observed pattern of neutrino masses and mixing. We show that, with a restricted 
number of such interaction terms with the A'-type couplings in the range (0.1-1.0), it is possible 
to forbid one-loop contributions to the neutrino mass matrix. This is illustrated with the help 
of a 'working example' where an econnomic choice of SUSY parameters is made, with three non- 
vanishing and 'large' R-parity violating terms in the superpotential. The two-loop contributions in 
such a case can not only generate the masses in the requisite order but can also lead us to specific 
allowed regions of the parameter space. 
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Neutrinos are massless to all orders in perturbation theory in the standard model (SM). However, the ever- 
accumulating data on solar, atmospheric and reactor neutrinos challenge us with the inescapable fact that neutrinos 
are massive and their physical states are mixtures of the flavour eigenstates [l|, [H, H[ . The SM has to be extended 
for explaining this. The simplest extension is the inclusion of 'sterile' right-handed neutrinos, whereby neutrinos 
may either acquire just Dirac masses or, with lepton number violation, participate in the see-saw mechanism which 
accounts for their ultra-light character. 

An alternative mechanism is provided by the supersymmetric (SUSY) extension of the SM with renormalizable R- 
parity (R p ) violating terms in the Lagrangian 0, @. The fact that baryon and lepton numbers are but accidentally 
conserved in the SM entails the possibility of R p = (— \^ B + L + 2S being violated in SUSY, where B, L and S are 
baryon number, lepton number and spin respectively. In order to avoid unacceptably fast proton decay, either B or L 
must be conserved, while the other may be violated. In the latter situation, small Majorana mass terms for neutrinos 
(with AL — 2) are generated, without the requirement of any additional fields Thus, the neutrino sector may be 
looked upon as a motivation for such L-violating interactions. 

The multiplicative conservation of i?-parity prevents the lightest SUSY particle (LSP) from decaying, as R p equals +1 
for all SM particles and —1 for the superparticles. All possibilities of i? p -violation are encapsulated in the following 
terms of the superpotential: 



where the first two trilinear terms and the bilinear term are AL — 1 and the third term is AB = 1. Since we are 
interested in neutrino masses, let us assume that B is conserved, and that R p is broken through L-violating couplings 
only. Moreover, we are neglecting the bilinear terms e^Li^ (on which we will comment later), and consider the 
trilinear A^ fc -type couplings only to illustrate our point. Elaborate studies in the recent years have led to constraints 
at various levels on these couplings [1, 0]- The pertinent gauge- invariant terms trilinear in particle/sparticle fields are 
given by 
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It is easy to see from above that the A^- fc -type couplings (27 of them altogether) can generate neutrino masses at the 
loop level, where the largest contribution comes from A^ 33 . We expect that all of the entries in the neutrino mass 
matrix should lie well within 1 eV. A generic expression for one-loop masses generated in this fashion is [8| 

K- lo °% ~ ^ m d k m d p M SU sY^Kk P ^pk, (3) 

where mf is the down-type quark mass of fc th generation, m 2 ~ is the (average) squark mass squared, and Msusy (~ Mj 
the Higgsino mass parameter) is the effective scale of SUSY breaking. If the masses thus induced have to answer to 
the observed pattern, then a SUSY breaking mass scale of about 500 GeV would in general imply A' ~ 10~ 5 — 10~ 4 
A similar conclusion follows for A-type terms, too. 

The question to ask is: are all trilinear R-parity violating couplings thus destined to be so small, irrespective of all 
other phenomenological considerations? For example, will the observation of any process which requires large values 
of some A'-terms mean that we need some additional mechanism to explain the neutrino mass pattern? We wish to 
demonstrate in this paper that it is not so, so long as one can eliminate the one-loop contributions but allow two-loop 
ones, through a limited number of A^-j, -terms. This drastically reduces the number of the A' terms whose signals may 
be of interest at the Large Hadron Collider (LHC), be it for direct observation or through indirect radiative effects. 

Situations where i? p -violating two-loop effects can contribute substantially compared to those at the one-loop level 
have been studied in earlier works 0]. In contrast, let us assume here a scenario in which there is a 'minimal' set of 
non-zero large (~ 0.1 — 1.0) A'- type couplings at the weak scale. One can clearly see from Eq. ^ that for such large 
A"s, it is impossible to explain the existing neutrino data, without going into unrealistically high values for m|, if both 
Kk P aim ^'j P k are allowed for the relevant {ij}-sets. A way out of this problem would be to postulate this minimal set 
of large A"s, of such composition that the above combinations do not exist, and the relevant interaction terms of Eq. 
(0 contribute to the neutrino mass matrix at the two-loop level (and beyond) only. At this level, together with the 
usual loop suppression factors with respect to the one-loop contributions, there will be additional sup pres sion coming 
from the parameters describing left-right mixing among different flavours in the squark mass matrices [lCj 1 . It is thus 
interesting to see whether these two suppression factors together may offset the 'largeness' of A"s, ultimately yielding 
contributions to neutrino mass matrix in accordance with the existing neutrino data. 

While there are no direct evidences of the Nature favouring any particular i?p-violating coupling over the others, one 
may, as a starting point, take those that are supported by the low-energy data. As a case in point, it has recently 
been advocated in [ll| that a minimal set of three i? p -violating couplings can simultaneously explain two interesting 
observations in flavour physics. The first one, as shown by the UTfit Collaboration, is the existence of a sizable 
deviation of the B s -B s mixing phase, j3 s , from its SM expectation, which is close to zero. The second one is the 
abnormally large leptonic branching ratio of the D s meson [l2j]. In ref. [ll|, it was found that one must have large 
A223 (A323) t° explain the recent data on D s — > fiv (D s — > tv) [13j], and in addition either \' 2 i2 or A 312 on similar 
order, contributing to the phase in B s mixing. The D s anomaly stems from a very accurate determination of the 
decay constant, Jd s , on the lattice by the HPQCD Collaboration 14j]- It has been pointed out in [l5[ that further 
clarifications are needed on some of the approximations used in [l4|, and prior to that, it may be advisable to use a 
more conservative estimate of Jd s , namely, (250± 15) MeV. Such a value is not in direct conflict with the experimental 
number (273 ± 10) MeV, and if one wishes to invoke i?j,-violation to explain the slight excess, one may use smaller 
values of the relevant couplings than those used in [llj. On the whole, we take the above result as a motivating 
feature of our analysis, without committing ourselves too decidedly on any specific numerical values. 

It may be in order to spell out at this stage how general our approach is, by re-iterating its main motivation. We 
would like to emphasize that we are not just attempting to compute two-loop diagrams contributing to neutrino 
masses, which have not been evaluated before. Nor is the sole purpose of this investigation to account for the claims 



1 Squark mixing parameters can in general occur along squark propagators, and may enter into one-loop contributions as well. We will 
talk about such loops in section [IT] In our analysis however, we have disallowed such combinations of A"s, and have only retained those 
which generate neutrino mass terms at the two-loop level. 
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on D s decays. The principal point made by us is that one can reconcile large R-parity violating couplings and neutrino 
masses, if only a subset of all possible couplings of such nature exist. If there is indication of large couplings, the subset 
must further be determined by the impossibility of generating one-loop neutrino masses. Two-loop contributions are 
tenable in such situations, and they can fit the entire neutrino mass matrix answering to the experimental constraints. 
An essential additional ingredient of this mechanism is SUSY flavour violation through squark mass matrices. We 
have stressed on identifying the minimum possible number of i? p -and flavour-violating parameters. This in a way 
restricts the set of contributing diagrams, but this feature is characteristic of a minimal choice and not of the specific 
couplings chosen, especially if the sfermions of different flavours are of comparable mass, a feature well-motivated 
from the suppression of flavour-changing processes. Thus this study reflects an entire set of possibilities rather than 
the property of some specific R-parity violating couplings. Let us also mention that the values of all R p - violating 
couplings are taken to be those at the electroweak scale and in the mass eigenbasis of the quarks. 

The paper is arranged as follows. In section [TT] we discuss the overall requirements in generating neutrino masses 
at the two-loop level only, using R p violating couplings of the A'-type. Some features of the two-loop contributions 
are outlined in section ITTTl while in section HVl we test the validity of the scheme on a numerical basis. Section fVl is 
on some correlated signals of the relevant couplings that may be tested at the LHC. We summarize and conclude in 
section PVTl Some representative expressions related to the loop integrals are included in the Appendix. 



II. THE PARAMETERS RELEVANT FOR TWO-LOOP EFFECTS 



Let us try to identify a minimal set of parameters that are required to generate a neutrino mass matrix at no less 
than the two-loop level. Of course, one requires a set of non-zero A' which can be allowed to lie in the range 0.1 — I. 
As will be explained below, one further requires the parameters controlling flavour violation in the squark sector in 
order to generate the mass matrix in a way consistent with observations. 

Next, we recall the pattern of the three-family neutrino mass matrix in the flavour basis, assuming, without any loss 
of generality, that the charged lepton mass matrix is diagonal in this basis. The constraints on the mixing angles are 

m 

sin 2 (26»i 2 ) = 0.86±g;gf 6 12 = (33.89 ± 1.44)° , sin 2 (20 23 ) > 0.92 =*> 9 23 > 36.8° , (4) 
and sin 2 (2#i 2 ) < 0.19. We assume the bilarge mixing scheme so that 823 = 7r/4 and #13 = [l7| . 



M v = 



I m 1 c 2 + m 2 s 2 ^§(-mi+m 2 ) - m 2 ) 

^(-mi + ma) \(m x s 2 + m 2 c 2 + m 3 ) \{-mis 2 - m 2 c 2 + m 3 ) ] , (5) 
^^(mi— m 2 ) \{-mis 2 - m 2 c 2 + m 3 ) \{mis 2 + m 2 c 2 + m 3 ) 



where mi, m 2 , m 3 are the mass eigenvalues, and s = sin #12, c = cos #12, 9ij being the mixing angle between the ith 
and the jth family. From this matrix one can easily take up the specific scenarios of normal (7773 >> 7772 > mi) or 
inverted (7772 > mi >> 7773) hierarchy or that of degenerate neutrinos (mi ~ 7772 ~ m 3 ). One can take mi = for 
normal hierarchy (NH) and 7773 = for inverted hierarchy (IH), without any loss of generality. In the case of NH, the 
existing data require (at 95% confidence limit) 



77; 



% = (7.60 ±0.35) x 10" 5 eV 2 , ImJ - m\\ = (2.50 ± 0.27) x 10~ 3 eV 2 , (6) 



and s 2 = 0.3. The corresponding numbers for IH and degenerate neutrinos (DN) are 

m\ - m\ = (7.60 ± 0.35) x 10~ 5 eV 2 , |m 2 - m\\ ~ \m\ - m§| = (2.50 ± 0.27) x 10" 3 eV 2 



(7) 



and 

mi ~ m 2 ^ m 3 ~ O(10 _1 ) eV (8) 
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FIG. 1: A typical one-loop diagram contributing to neutrino masses. It should be remembered that corresponding to each 
such diagram there is one with Vi and Vj flipped. 



respectively [l8| . 

Let us first try to understand intuitively the properties of the 'minimal set'. 

• There must be no less than three A' type couplings, each with a different leptonic index, for the three neutrinos. 

• To prevent mass generation at one- loop, couplings like A^ , which generate diagonal entries of the neutrino mass 
matrix, are forbidden. 

• Similarly, combinations like A^ fc; A^ fc are forbidden to prevent the off-diagonal entries at the one-loop level. 

• In fact, X' ikl \j mk combinations are also not allowed, since they can generate one-loop masses with the mass 
insertion S^f. 

This leaves us with a limited number of possible choices. 

As already discussed in Section HI our choice of the supersymmetric scenario is partially motivated by the explanation 
of the results on D s decays 2 . We thus include A 223 and A 323 in our minimal set of A'-type couplings, and propose 
that their values be allowed to be large, consistent with the individual constraints. It is easy to see from the relevant 
interaction terms (the second and third terms of Eq. |(5J)) that we need one more A£- fe with i = 1, in order to have 
contributions to the elements in the first row and the first column of M v . The choices that we thus have are A' 112 , 
^m> ^n3) -^131 an d A' 123 . 

Let us clarify the last criterion mentioned above. As a first choice, let us choose A' 112 . It is then easy to see from Figure 
1 that there are non-vanishing one-loop contributions to the (1,2) and (1,3) elements of the neutrino mass matrix. 
This is because of the fact that the quark and squark mass matrices of the same charge are not in general diagonal 
simultaneously; the evolution of the squark mass parameters from the high scale of SUSY breaking always tend to 
destroy such alignment. The resulting possibility of a flavour transition as well as a chirality flip along the down- type 
squark propagator allows one to obtain some one-loop contributions to the neutrino mass matrix. These contributions 
are driven by a parameter <5f 3 fi of mass-squared dimension, which is basically the corresponding off-diagonal term in 
the down-type squark mass matrix 3 . Such diagrams are not suppressed enough to balance the large values (O(0. 1)) 
of A 223 and A 323 and give admissibly small entries for the (1,2) and (1,3) elements of the neutrino mass matrix. 
So, a non-zero A' 112 will not normally serve our purpose. Besides, the phenomenological constraint on A' 112 makes it 
an inappropriate candidate for the demonstration of the effects of large R p violating interactions. Following similar 
arguments, the choice of A' 131 should be abandoned in our minimal set of A"s. 

On the contrary, none among A' 121 , A' 113 and A' 123 can give one-loop contributions to the neutrino mass matrix. So 
in principle any one of them can be included in the minimal set together with A 223 and A 323 , to generate neutrino 



2 This is a partial motivation because, as we will show later, the allowed values of Aj 2 3 and A 323 result only in a marginal enhancement 
of the D s leptonic branching ratio. However, it is better to be cautious about the HPQCD lattice result. 

li Our convention is different from, say, that of While our 8, which is of mass-squared dimension, is identical to their A, the A 

parameters that we subsequently introduce are based on a different scaling. We have checked that the existing numerical constraints 
are all satisfied. 
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mass at two-loop level. Note that the choice of A' 123 puts a single squark mixing parameter at our disposal, namely, 
the one describing the second-and third- family squark mixing (#23 )■ Thus we would have a set of four independent 
parameters: A' 22 3, A 323 , A' 123 and S^P". However, as will be evident from our numerical results in Section llVl it is 
difficult to fit the six independent elements of M v with experimental data with just these four parameters. 

Choosing A' 113 , on the other hand, will involve ^f 3 fl , the first- and third-family squark mixing parameter, for generating 
the elements of M v in the first row and first column, and S23 for the rest of the matrix elements (from now on, we will 
drop the chirality superscript on the 6s, since the only type that we will ever be interested in are those of the LR type 
in the down-squark sector). This means that for this choice we have a set of five independent parameters comprising 
°f (^223> ^323> ^ii3> ^23 and #13). In Section HVl we will see that in this case we are able to fit elements of M„ with the 
existing constraints. In a similar way, the choice of A' 121 also leads to the same number of independent parameters. 
However, for the latter choice, some two-loop contributions would be suppressed further by the ratio m s /nib, making 
the two- loop effects undesirably small, as we shall see in Section [TV] 

Thus, our selected parameter space consists of a minimal set of three O(0. 1) A"s, namely A' 223 , A 323 and A' 113 , and 
two non-zero squark mixing parameters 613 and 623, generating neutrino masses at the two- loop level. All other 
parameters are set to be zero at the weak scale. Also, we will work under the assumption of all the A"s being real. 

In our calculation, we scale the squark mixing parameter Sij by the factor m^MsusY , and define a dimensionless 
parameter A^- = &y / '{m^M su sy) , with the already specified connotation Ay — A^. The various loop contributions 
which involve flavour violation and require a chirality flip in the (down-type) squark propagator are expressed in terms 
of Ai 3 and A 23 . 

The coupling A' 113 is bounded from charged current universality [l9| as well as processes like n + — > e + v e . Here we 
use a 99% confidence level bound of |A' 113 < 0.15. As we shall see, this relative smallness of A' 113 leads to a distinct 
preference of the NH scenario of neutrino masses over IH or DN. 

A 223 and A 323 , the other two couplings, can be large, even 0(1). We have checked that the recent CLEO constraint 
on lepton flavour violation in T — > /it [20| is consistent with this upper limit. 




Ref. 2 of [6( quotes a weak scale bounds of |A'n 3 |, A'i 23 | < 0.39. These limits arise from the need to prevent tachyonic 
sneutrinos even at the Grand Unified Theory (GUT) scale [211 ] . The maximum value at the GUT scale is driven by 
the input parameters; for the set known as SPSla, this comes out to be about 0.13. When run down at the Mz scale, 
the coupling increases threefold and the bound becomes 0.39. One can easily relax this bound for other choices of the 
GUT scale input parameters. 

The dimensionless parameters Ajj can be constrained from various flavour-changing neutral current (FCNC) processes. 
For those that we are interested in, Ai 3 is constrained from B°-B° mixing to be less than 5.2, and A 23 is constrained 
from the inclusive b — > sj branching ratio to be less than 1.0. 

Let us mention again that this is just one of several possible choices. Following the rules laid down earlier, one must 
have three A' type couplings and two <5-type squark mixing parameters. However, some of the possible choices are 
extremely constrained from data. For example, the choice of A' 121 , A 221 , A 323 , 6^ and 8^ is severely restricted by 
the absence of leptonic flavour-violating decays ir° — > e/i, <fi — > e[i, B — > e(/x)r etc. 



Having shown that there are no one- loop contributions to the neutrino mass matrix M v , let us enlist and compute 
the two-loop contributions that are driven by the three nonzero A' type couplings and two 8 parameters. We work in 
the 't Hooft-Feynman gauge. 

Although the individual diagrams are divergent, the very fact that there is no counterterm at the tree-level for the 
interactions generated at higher loop levels immediately tells us that the end result is finite. This is ensured when all 
diagrams including all possible fields and their superpartners are taken into account. 




III. THE TWO-LOOP CONTRIBUTIONS 
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FIG. 2: Two-loop diagrams that make leading contributions to neutrino masses. The contributions in the three different 
diagrams are proportional to (a) X' ipk X' jk , p , A^V Un d h V* ndk , (b) X' ipk X' jk , p , A^V Uk , dk (c) X' ipk X' jk , p , A^K p d fe , ■ For our choice 
of A''s in (c), k — p' = 3 and thus A kp i = 1, so that there is no squark flavour violation. The flipped diagrams, too, will 
contribute as usual. 



Figure 2 represents three classes of diagrams that turn out to be dominant in our study. A full list of generic expressions 
for the loop-factors arising from these diagrams is provided in the Appendix. The amplitude corresponding to diagram 
2(a) for the (1, 1) and (i, j) elements of M v , where i, j = 2,3, is found to be 

respectively, where £ t = V t * s V t b, and £' t = V t * d V t b- The loop functions have been left out of these expressions. For the 
(1, 2) and (1, 3) elements of M„, the contribution from diagram 2(a) contains two separate parts proportional to the 
two factors written above, along with the appropriate loop functions multiplying each of them. 

The amplitude for diagram 2(b) vanishes when mediated by W ± (but not the charged Higgs or Goldstonc field), 
which follows from the details of 7-matrix algebra. For the diagonal entries of the neutrino mass matrix, this diagram 
yields 



Mil ~ V ub ^A 13 x e x u , 
M w 

lvl w 

K 2b ] 33 - V cb ^A 23 x T x c (10) 



where x a — m 2 a /M^ v for a — e, yu, r, u, c. Each of the off-diagonal matrix elements (1,2), (1,3) and (2,3) is a sum of two 
terms which are respectively proportional to the first and second, first and third and second and third factors written 
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above, again with the corresponding loop functions. In general, being proportional to the squares of lepton masses, the 
contribution of diagram 2(b) are suppressed compared to those of Figure 2(a). Diagram 2(c) is particularly interesting, 
since there is no flavour change required along the internal squark-line in this diagram. Thus, these diagrams do not 
have the suppression by A-factors. Nevertheless, this diagram has an overall lepton mass dependence. So ultimately 
it contributes more than diagram 2(b), but less than 2(a). For the diagonal entries of the neutrino mass matrix from 
2(c), 

[M^lu ~ V; d x u m e , 
[M 2 V C } 22 ~ V*x c m^, 

[ M " C ] 33 ~ V*x c m T (11) 

while, just as before, each of the off-diagonal entries M u (i,j) separately contains two terms which are proportional 
to the i-th and j-th factors respectively of Eq. (TTT]) . It is thus clear from equations ©, (flU)) , and (fTTj) , that the 
contribution from diagrams of type 2(a) dominate over the others for the elements in the first row and first column 
of M v , while for the other elements, these are more or less of the same order. 

We have worked with such a choice of the electroweak symmetry breaking sector that tan (3 = 10 (where tan (3 is 
the ratio of the two Higgs vacuum expectation values) and the charged Higgs mass is 500 GeV. The charged Higgs 
contributions are found to be suppressed with respect to the ones discussed above. In a similar manner, the loops 
involving charginos and neutralinos are found to be of subleading nature, as their presence would imply additional 
squark and slepton propagators, leading to bigger suppression factors under our choice of mass (~500 GeV) for all 
squarks and sleptons. It is therefore legitimate to illustrate our main points leaving out such diagrams. Additional 
diagrams have been taken into account in earlier works dealing with two-loop neutrino masses in R-parity violating 
SUSY Hj. Representative diagrams of this type are shown in Figure 3. The reasons for not taking these contributions 
into account, without losing generality in our approach, are as follows: 



• The contributions from diagrams 3(a) depend on the splitting between CP-even and CP-odd sneutrino states. 
That requires added theoretical inputs which are not present in our study. 

• Even when one goes beyond the minimal set of R-parity violating interactions, diagrams of the kind shown in 
3(e) cannot contribute without there being contributions at one-loop, whose absence is precisely the theme of 
our work. 

• Diagrams 3(f) require additional assumptions about soft trilinear terms with AL = 1 in the scalar potential. 



IV. RESULTS AND DISCUSSION 



We have five parameters, namely, A' 113 , A 223 , A 323 , A23 and A13, with which to fit the neutrino mass matrix M„ to 
generate the required mass hierarchies. Here, as we have already defined, Ay = 8ij /m^MsusY ■ These parameters, 
along with 

m t = 172.5 GeV, m 6 = 4.5GeV, \V U \ = (8.12 ± 0.88) x 10~ 3 , \V ts \ = (40.67 ± 1.30) x 10~ 3 , 

\V cb \ = (40.8 ±0.6) x 10~ 3 , sin 2/3 d = 0.755 ±0.040, 12 = (33.89 ± 1.44)° , (12) 

(where (3d — arg(V^)), all sparticle masses (including H^, and all sleptons and squarks, and hence M$usy) at 500 
GeV, and tan (3 = 10, essentially determine the entries of M v as shown in the Appendix. 

We vary the SM inputs over their allowed ranges, and the five parameters A' and A over the range 0.0-1.0, to see 
whether any simultaneous solution to the M„ constraints exist. We take all the A's and As to be real and positive. In 
fact, there are four independent parameters, and not five, that need to be varied. The reason lies in the neutrino mass 
matrix Af„, whose (2,2) are (3,3) elements are equal for #13 = and differ very slightly for small f?i 3 . The relevant 




FIG. 3: Additional two- loop diagram that will not contribute in our case, (b), (c), (d), correspond to the blob shown in (a). 
Contributions from (/) require trilinear //-violating soft terms in the scalar potential. 



amplitudes, being completely identical in the leptonic part, imply A' 22 3 ~ A 323 . Thus, essentially, we have four free 
parameters, namely, A' 113 , A 223 , A13 and A 2 3. 

As a result of varying all the parameters, there are six possible projections of the four-dimensional scatter plot. In 
figures HE we show four of them, the other two not giving any independent information. 

The plots are drawn for (i) NH (fig-H, (ii) IH (fig. [5), (iii) DN (fig.©, all with 6 13 = 0, and (iv) NH (fig. EJ), with 
#13 = 10°. No such figures are separately shown for the IH and the DN cases, because (a) there is no appreciable 
difference with the corresponding #13 = case, and (b) these scenarios are in general disfavoured by the constraints 
on A' U3 . 

The scatter plots are essentially based on the fact that, corresponding to any value of one of the four aforementioned 
parameters, we get confined to rather narrow intervals of the remaining three, in order to satisfy the relative values 
of the neutrino mass matrix elements answering to the tri-bimaximal mixing pattern. Thus the scatter plots turn 
into correlation curves (bands) whose widths come largely from the uncertainties of the neutrino oscillation data and 
marginally from the uncertainties in the CKM elements. 
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FIG. 4: Correlation plots for NH case with 8ri = 0. The vertical (horizontal) line in the top right (bottom left) panel corresponds 
to the 99% CL upper limit on X' 113 (see text). 



To see the most important conclusions, let us first concentrate on fig. 2] The upper panels show the allowed regions 
for A' versus A; one goes up as the other goes down. Qualitatively, this can be understood from the expressions of the 
M„ elements as given in the appendix: the product of the type A'A'A appears in the leading contributions. The lower 
left-hand panel shows the correlation between A' 113 and A 223 ; taken in conjunction with the upper panels, this also 
tells the allowed regions of the corresponding As, and this fact has been confirmed in the lower right-hand panel. The 
upper bound on A' 113 , shown by a vertical line in the upper right-hand and by a horizontal line in the lower left-hand 
panels, corresponds to the 99% confidence level limit from charged current universality. 

It should be noted that the allowed regions for A' 113 fall outside this limit for both IH and DN cases. However, 
such constraint, as listed in existing literature, assumes the existence of no A-type couplings, which can invalidate 
the bound but play an ineffective role in neutrino mass generation, giving contributions suppressed by light lepton 
masses. With these in view, we have allowed A' 113 to have values larger than the upper bound found in the literature, 
with the caveat that the large values may indicate the existence of additional interactions of A-type. 

The graphs clearly show that the NH scenario favours larger values of A 223 than IH, while for A' 113 it is the other way 
around. This is because, in the IH case, one requires the (1, 1) element of M v to be of higher magnitude, and one is 
at a relative disadvantage in the loop contributions, since the contribution to this element is suppressed by the down 
quark mass. One also gets restricted to rather small values of A 223 is this case. For the degenerate neutrino case, too, 
A' 113 has to be on the higher side, since the corresponding contributions do not get the advantage of heavier quark 
masses. This re-iterates the difficulty in reconciling the IH and DN scenarios with the constraints on A' 113 , which can 
be bypassed through, for example, the occurrence of additional R p violating interactions. 
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While it is true that the preference of NH over IH results from the way we have selected our parameters, it should be 
also be noted that it is more the result of selecting three X' -type couplings in the range of 0.1 and two squark flavour- 
violating parameters A. It is of course true that one can fit the IH and DN scenarios with a larger set of R-parity 
violating interactions. However, with the so-called 'minimal' choice, the orders of magnitudes of the loop contributions 
are not significantly different, so long as the A'-parameters are in the same range, and the squarks of different flavours 
(due to the different indices of these parameters) participating in the loops are in the same mass range. Naturally, the 
contributions will be smaller with more massive squarks; then higher values of the R-parity violating couplings than 
what is indicated can be accommodated. Thus, while it is not our goal to establish the preference of one scenario over 
the other, what we successfully show is that one can generate neutrino masses with large R-parity violating couplings, 
and that a pattern follows from a minimal choice, which does not necessarily depend on which three parameters are 
involved. 



V. CORRELATED SIGNALS: SOME SPECULATIONS 



As we have noted earlier, the small values of A' 22 3 and A 323 can enhance the D s — > [i(j)v branching ratio marginally. 
However, if one indeed entertains the possibility of some other A-type interaction to save the IH or the DN picture, it is 
possible that these two couplings may become large. The lepton flavour violating (LFV) decay T — > /jt is, again, only 
marginally enhanced, and is still well below the experimental limit. However, a positive signal in this channel would 
be very interesting from the neutrino perspective. The same comment applies to other LFV decays, like D° — > e/i, 
driven by A' 113 A 223 . 
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One of the most interesting low-energy effects for this scenario is the change in the branching ratio of K + — > ■n+vv. 
The decay, based on s — > dvP, is again controlled by A' 113 A2 2 3- The experimental number is B(K + — ► tt + i>D) = 
(1.47±o;|g) x 10~ 10 16], while the SM prediction is about (0.8 ±0.1) x 10" 10 . It was pointed out in 22] that an exact 
upper bound is difficult to obtain considering the interplay of the SM, the R p conserving SUSY and the R p violating 
SUSY, but it can safely be said that with couplings of the order that we have used in this work, the R p violating 
amplitude may even be larger than the SM amplitude. In that case, this mode cannot be used as a clean channel for 
extracting sin(2/3). Measurement of the said angle and a comparison with the charmonium result will again be crucial 
for our ansatz. 



VI. SUMMARY AND CONCLUSIONS 



We have considered scenarios where R p violating couplings can be large, and the neutrino mass matrix can still be 
generated in a manner consistent with observed results. This, we argue, can be possible if there are only a few 
couplings of this type, so that the combinations necessary for one-loop neutrino masses are not available. Two-loop 
contributions come to one's advantage in such situations, together with the possibility of flavour violation in the 
sfcrmion mass matrices. Considering the A'-type couplings, we have demonstrated this; with three such couplings 
and two squark flavour violating parameters, the NH scenario can be reproduced, guiding one to a specific region of 
the parameter space. For the IH and DN cases, however, this requires the value of at least one coupling to come into 
conflict with observable constraints unless one postulates additional R-parity violating terms in the superpotential. 
Of course, there may be more than one choice of the set of R-parity violating couplings leading to two-loop neutrino 
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masses, and the exact numerical consequences in the neutrino sector can be dependent on which A'-couplings actually 
exist. 

Similar conclusions can be established if one includes the bilinear R-parity violating terms in the superpotential. One 
neutrino state acquires a tree-level mass in such a case, thus relaxing the constraint that seems to loom large on the 
parameter A' 113 as discussed above. The two remaining couplings (with values in the range 0.1-1.0) and the squark 
flavour violation parameters can then generate the remaining terms in the mass matrix at the two-loop level. This 
may make the IH and DN cases less constrained. 

In conclusion, large trilincar R-parity violating interactions are not necessarily an impediment to the explanation of 
neutrino masses and mixing. Thus if some phenomenon observed in the laboratory points strongly towards such large 
interaction strength, it may still explain the neutrino sector perfectly well, provided that only a few R-parity violating 
interactions occur in nature with sizable strength. 
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The two- loop matrix elements are expressed in terms of the following variables. 



6 - KVn 
3 = v; d v tb 
£ c = v; s v cb 
£ = v; d v cb 

x t = mf/M^ 



(13) 
(14) 
(15) 
(16) 
(17) 



The generic loop-functions with proper arguments are listed below. There are two types of such functions, depending 
on whether they depend on lepton masses or not. 



Functions, first set: i = e,fj,,T. 



F^xuxw) = 



3x t -l x t 2 logx t 



4(x t - 1) 2{x t - l) 2 



(18) 



F 2 (x t ,x w ) = 1 



xt log St 
x t -l 



(19) 



Fi(x t ,x w ) - F 2 (x t ,x w ) = ^} Xt \ogx t 



2(x t -l) 2 b 4(^-1) 



(20) 



Fi(x c ,Xi) 



F 3 (x t ,x w ) = -2Fi(x t ,xw) 



1 - x c + x c log Xc_ Xc\ogx c — x c -Xi log Xi + Xi 



1 - x c 



X c Xi 



(21) 
(22) 



^5 (^c? %i) 



1 rl 2 _ l 2 _ 



(x c - x w ) 2 

1 2 , 3 

--x w logx^ - - 



(23) 



l 



(x w 3?c) 



ix^ log x w - ^x 2 w - ix 2 log x c - Jx 2 

-(x w <-> Xi) 



(24) 



F 7 (x c ,Xi) = 



l r,l 



■( 



1 3 ^ 2 2 \ 1 u , "^"^3 

g-^M) _ ~^ X W X C T X W X C ) — q X C 1°§> X C ~ Yg X ° 



(X W Xi) 



Functions, second set: 



F 8 (Xq,X b ) 



x 6 (logx^ - logx b ) 



(Xq - X b ) 2 



(25) 



(26) 



F 9 (Xq,X b ) = 



(Xq - X b ) 2 



ix? log Xq-^x\- X q X b log Xq + X q X b + \^x\ log X b - ^x\ 



(27) 



F 1Q (Xg,X b ,X s ) = F g (Xg,X b ) - F g (Xg,X s ) 
F U (Xq,X b ,X s ) = F 6 (Xq,X b ) - F 6 (Xq,X s ) 



Fl 2 {Xq,X b ) = 



log X q - log X b 

Xq Xf) 



F 13 (Xg,X b ) = 



(X q - X b ) 



\x\ logxg - \x\ - \x\ logxt + \x\ 



Matrix element M v (l, 1): 
M v (l, 1): Diagram 2(a) with </> 



Aii 3 Aii 3 4M ^ (167r2)2 



C> M mbmd A 
O M SUS y — 2 2 A 13 

^5 — 



m 2 (Fi - i^Fio^,^,^) 



-m 2 t x d F 2 Fn{xq, x b , x d ) + m d F 1 F 10 (xg, x b , x d ) 



M v (l,l): Diagram 2(a) with W 



, g 1 m b md r 

A 113 A 113 T 7TSZ2\2 Msus V ^2 — A 13 



4 (16^ 2 ) 2 to 2 - m 2 



£t [F3 (xt , x w ) - F 3 (x u , x w )] 



+C[ F 3 (x c ,X W ) -F 3 (x u ,x w )\ F W (Xq,X b ,X d ) 
M„(l, 1): Diagram 2(b) with <p 



X' X' 

A 113 A 113 



4M 2 , (16^ 2 ) 2 



m b m d A13 



TYl~ 

9 9 [Xu-^4 \Xu , Xe)t*% \Xq, X b ) 

w - m i 



-F 4 (x u ,X e )F 9 (Xq,X b ) - F 5 (x u ,X e )F 9 (Xq,X b )] 

M v (l, 1): Diagram 2(c) with </> 



A 113 A 113 Me-^ V ud 



(167T 2 ) 5 



™ e F 5 (x u ,a; e )Fi3(x ( j,a; (i ) 



.X, 



■m e [F 7 (x u ,x e )F- L3 (xg,x d ) - F 6 (or u ,a; e )i ? i2(a; 5 ,a; ( i)] 



Matrix clement M„(2,2): 
M„(2,2): Diagram 2(a) with 



X' X 



.9 



1 



m b m s 



6 M, 

mfx s F 2 F 11 (xq, Xb, x s ) + m 2 s FiF w {x q , x b , x s ) 



-2,-223 ^ (167r2)2 M - raaB m2 _ m2 ^23 



m 2 (Fi - F 2 )Fi (a;g, ir 6 , x s ) 



M„(2,2): Diagram 2(a) with 



X' X' 
A 22 3 A 22 3 



4 (16tt 2 ) 



m b m s 



2\2 ±VJ -su8 V 2 



mf - mi 



i 2 :5 



6 [-F3 (at , «w ) - ^3 (a;« , % )] 



+£c[F 3 {x c ,X W ) - F 3 (X U ,X W )) F W (Xq,Xb,X s ) 



M„(2,2): Diagram 2(b) with 



X' X' 

A99Q Ac 



223 A 223 



4M 2 , (167T 2 ) 2 



V c b M sus ,, m b m s A 23 



MS 



-F4(.T c ,a; M )F9(a;q,x h ) - F 5 (x c , x^)F 9 (x q , x b )} 
M v (2,2): Diagram 2(c) with 



X' X' n 2 

A 223 A 223 y 



(16tt 2 ) 



2\2 cs 



+- 



X w X , 



-mn[F 7 (xc,Xn)Fi 3 (xq,Xs) - F 6 (x c , x^F^xg, x s )] 



Matrix element M u (3, 3) same as M v {2, 2) with fi replaced by r. 



Matrix clement M„(l, 2) = M„(2, 1): 
M 1/ (l,2): Diagram 2(a) with 



X' X' 

A 113 A 223 ofl^2 /ic„2\2 -"*s«siM 



8M 2 , (16tT 2 ) 5 



M s . 



« ^2 A 13 



m 2 (Fi - F 2 )Fi (xq, ir s ) 



-m 2 x s F 2 Fu(x q , x b , x s ) + m^FxFwixq, x h , x s ) 



+ 



, m b m s - 
St — 2 2 A23 



m 2 (Fi - F 2 )Fio(a;g,a;6,2;d) 



-m 2 x d F 2 F 11 (xq,x b ,x d ) +m 2 ,F 1 F w (xq,x b ,x d ) ) 
M„(l,2): Diagram 2(a) with 



A 'll3 A 223 ~ TTaZrVl Musi/ 



m b m d 



— 2 2 ^13 

m„ — mf 



£t[F 3 {x t ,x w ) - F 3 (x u ,x w )] 



m b m s 



2 A 23 



8 (16tt 2 ) 

+£c[F 3 {x c , x w ) - F 3 (x u ,x w )] F 10 (xq, x b , x s ) 

-F 3 {x u ,x w )\ +£,' c [F 3 {x c ,x w ) - F 3 (x u ,x w )] F w (xq,Xb,x d f) 
M,(l,2): Diagram 2(b) with <p 

9 2 1 



£' t [F 3 (x t ,x w ) 



Ai 1 o Xn 



M, 



113 223 8M W (167T 2 ) 2 ^ susy < 



V c6 m 6 m d Ai 3 



X c F4(x c ,X^)F 8 (Xq,X b ) 



-F i (x c ,x ft )F 9 (xg,x b ) - F 5 (a; c ,a: M )F9(a;5,a; fc )] 



+ 



W ""n 

V ub m b m 8 A 23 



—2—2 — 2 -[a;„i ; 4(a; u ,a; e )F 8 (x ( j,a;{ ) ) - F 4 (x u , x e )F 9 (x q , x b ) - F 5 (x„, x e )F 9 (xq, x b )] 
M w m e 



M„(l,2): Diagram 2(c) with 



X' X' a 2 
A 113 A 223 Si 



+- 



X \ 

-f m M [F 7 (a; c ,a; A1 )Fi3(x I j,a; s ) - F 6 (x c ,x^)F 12 (xq,x s )}) 

X 11 ' 



+ 



V. 



■ad 



•E'u %d 



+- 



X w X, 

x d 



m e F b (x u , x e )F 13 (xg, x d ) 



Xyj X( 



m e [F 7 (x u , x e )F 13 (xq, x d ) - F 6 (x u ,x e )F 12 (Xq,x d )fj 
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Matrix elements M v (l, 3) and M v (3, 1) are same as M„(l, 2) with /x replaced by r. 
Matrix element M„(2, 3) = M„(3, 2): 



M I/ (2,3): Diagram 2(a) with 



Ai)QQ Ao' 



6 M„ 



m b m s 



A 2 



susy 9 9 ^—^23 

771, — mr 



223 323 8M 2 , (16tt 2 ) 2 
-m^x s F 2 F 11 (xq, x b ,x s )+ m 2 s FiF 1Q (xq, x 6 , x s ) 



m 2 (Fi - F 2 )Fi (xq, x b , x s ) 



M„(2,3): Diagram 2(a) with W 



/\99q/\qoq — — 



m b m s - 

jyj-susy n 9 ^23 



^223-323 g (16?r2)2 m 2 _ ^ 

+^c[-F3(»c,a:w) ^ F 3 (x u ,x w )] F 10 (Xg,x b ,x s ) 
M„(2,3): Diagram 2(b) with <j> 



£t [Fa (x t , x w ) - F 3 (x u , x w )] 



A223A323 m ^ (167r2)2 



V c b M SUS y m b m s A 2 3 



- ml 



-F^XcX^FgiXq^Xh) - F 5 (x c , X fl )Fg(Xg, x b )] + 
-F 4l (x c ,X T )Fg(Xq,X b ) - F% (x c , X T )Fg [Xq, X b )\ 

M„(2,3): Diagram 2(c) with 



M^ w -ml 



[x c F 4 (x c , X^FsiXq, X b ) 



[XcF^Xc, X T )F 8 (Xq, X b ) 



\l \l 2 1 
^223^323 9 ( 167r 2)2 "cfc 



-m / ^F 5 (a; c ,a; (ti )-Fi3(a;g,a; s ) 



-m^[F7(x c , X,j,)Fl3(Xq, x s ) - F 6 (x c , X^)F 12 {Xq, x s )] 



-( c _ s m T F 5 (x c ,x T )F 13 (xq,x s ) 

\ Xw Xr 



H m T [F7(x c ,a; T )Fi3(x^,a; s ) - F 6 (x c , x T )F 12 (xq, x s )] 



Xt qjj 00 



(44) 



(45) 



(46) 



(47) 



[1] Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 87, 071301 (20 01) |arXiv:nucl-ex/0106015j ; Q. R. Ahmad et al. 
[SNO Collaboration], Phys. Rev. Lett. 89, 01130 1 (2002) [arXiv:nucl-ex/0204008| ; Q. R. Ahmad et al. [SNO Collaboration], 
Phys. Rev. Lett. 89, 011302 (2002) jarXiv:micl-ex/0204009] ~ 

[2] A. Bandyopadhyay, S. Choubey, S. Goswami and K. Kar, Phys. Lett. B 5 19, 83 (2001) arXiv:hep-ph/0106264j ; 
M. V. Garzelli and C. Giunti, JHEP 0112, 017 (2001) [arXiv:hep-ph/0108191 j ; J. N. Bahcall, M. C. Gonzalez-Garcia 
and C. Pena-Garay, JHEP 0108, 014 (2001 ) [arXiv:hep-ph /0 106258]; G. L. Fogli, E. Lisi, D. Montanino and A. Palazzo, 
Phys. Rev. D 6 4, 093007 (2001) j arXiv:hep-ph/ 0106247 ; P. L Krastev and A. Y. Smirnov, Phys. Rev. D 65, 073022 
(2002) [arXiv:hep-ph/0108177 ; A. Bandyopadhyay, S. Choubey, S. Goswami and D. P. Roy, Phys. Lett. B 540, 14 (2002) 
[arXiv:hep -ph/0204286 . 

[3] K. Eguchi et al. [KamLAND Collaboration], Phys. Rev. Lett. 90, 021802 (2003) [arXiv:hep-ex/0212021j . 

[4] G. R. Farrar and P. Fayet, Phys. Lett. B 76, 575 (1978); C. S. Aulakh and R. N. Mohapatra, Phys. Lett. B 119, 136 
(1982); L. J. Hall and M. Suzuki, Nucl. Phys. B 231, 419 (1984); J. R. Ellis, G. Gelmini, C. Jarlskog, G. G. Ross and 
J. W. F. Valle, Phys. Lett. B 150, 142 (1985); G. G. Ross and J. W. F. Valle, Phys. Lett. B 151, 375 (1985); S. Dawson, 
Nucl. Phys. B 261, 297 (1985); R. Barbieri and A. Masiero, Nucl. Phys. B 267, 679 (1986). 



17 



[5] S. Dimopoulos and L. J. Hall, Phys. Lett. B 196 (1987) 135; R. M. Godbol e, P. Roy and X. Tata, Nucl. Phys. B401, 

(1993) 67. R. M. Godbole, P. Roy and X. Tata, Nucl. Phys. B 401, 67 (1993) [arXiv:hep- ph/920925l] 
[6] B. C. Allanach, A. Dedes and H. K. Dreiner, Phys. Rev. D 60, 056002 (1999) arXiv:hep-ph/9902251 ; Phys. Rev. D 69, 

115002 (2004) [Erratum-ibid. D 72, 079902 (20 05)] [arXiv:hep-ph/0309196| . 
[7] See, for example, H. K. Dreiner, arXiv:hep-ph/9707435 G. Bhattacharyya, arXiv:hep-ph/9709395 R. Barbier et 

al, |arXiv:hep-ph/9810232] S. Rakshit, G. Bhattachar^a and A . Raychaudhuri, Phys. Rev. D 59, 091701 (1999) 

[arXiv:he p-ph/9811500 ; B. Mukhopadhyaya, |arXiv:hep-ph/0301278| 
[8] M. Drees, S. Pakvasa, X. Tata and T. ter Veldhuis, Phys. Rev. D 57, 5335 (1998) [arXiv:hep-ph/9712392| ; E. J. Chun, 

S. K. Kang, C. W. Kim and U. W. Lee, Nucl. Phys. B 544, 8 9 (1999) [arXiv :hep-ph/9807327 . 
[9] F. Borzumati and J. S. Lee, Phys. Rev. D 66, 115012 (2002) [arXiv:hep- ph/0207184|7^ 
[10] J. S. Hagelin, S. Kelley and T. Tanaka, Nucl. Phys. B 415, 293 (1994); S. Bertolini, F. Borzumati, A. Masiero and 

G. Ridolfi, Nucl. Phys. B 353, 591 (1991); F. Gabbiani, E. Gabrielli, A. Masiero and L. Silvestrini, Nucl. Phys. B 477, 

321 (1996) [arXivihep ph/9604387] ■ 

[11] A. Kundu and S. Nandi, Phys. Rev. D 78, 015009 (2008) [arXiv:0803. 18981 [hep-ph]]. 

[12] B. A. Dobrescu and A. S. Kronfeld, Phys. Rev. Lett. 100, 241802 (2008) [arXiv: 0803. 05121 [hep-ph]]. 

[13] C. Amsler et al. [Particle Da ta Group], Phys. Lett . B 667, 1 (2008); B. Aubert et al. [BABAR Collaboration], Phys. 
Rev. Lett. 98, 141801 (2007) [arXiv:hep-ex/0607094l ; T. K. Pedlar et al. [CLEO Collaboration], Phys. Rev. D 76, 072002 
(2007) [arXiv:0704.0437l [hep-exjj; K. M Ecklund et al. [CLEO Collaboration], Phys. Rev. Lett. 100, 161801 (2008) 
|arXiv:0712.TT75l [hep-ex]]; K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 100, 241801 (2008) |arXiv:0709.T340l 
[hep-ex]]. 

[14] E. Follana, C. T. H. Davies, G. P. Lepage and J. Shigemitsu [HPQCD Collaboration and UKQCD Collaboration], Phys. 

Rev. Lett. 100, 062002 (2008) arXiv:0706T726] [hep-lat]]. 
[15] V. Lubicz and C. Tarantino. rarXiv:0 807.4605 [hep-lat]; 
[16] See Ref. 1 of [Hj]. 

[17] P. F. Harrison, D. H. Perkins and W. G. Scott, Phys. Lett. B 530, 167 (2002) [arXiv:hep-ph/0202074| . 
[18] A. Bandyopadhyay, S. Choubey, S. Goswami, S. T. Petcov and D. P. Roy. larXiv:0804.4857l [hep-ph] . 
[19] V. D. Barger, G. F. Giudice and T. H an, Phys. Rev. D 40, 2987 (1989). 

[20] W. Love et al. [CLEO Collaboration], larXiv:0807.2695l [hep-ex], 

[21] B. de Car los and P. L. White, Ph ys. Rev. D 55, 4222 (1997) [arXiv:hep-ph /9609443| . 
[22] J. Welzel, |arXiv:hep-ph/0505094| 

[23] S. Roy and B. Mukhopadhyaya, Phys. Rev. D 55, 7020 (1997) [arXiv:hep-ph/9612447l ; B. Mukhopadhyaya, S. Roy 
and F. Vissani, Phys. Lett. B 443, 191 (1998) [arXiv:hep-ph/9808265| ; M. Hirsch, M. A. Dia z, W. Porod, J. C. Ro - 
mao and J. W. F. Valle, Phys. Rev. D 62, 113008 (2000) [Erratum-ibid. D 65, 119901 (2002)] [arXiv:hep-ph/0004115j ; 
A. S. Joshipura, R. D. Vaidya and S. K. Vempati, Nucl. Phys. B 639, 290 (20 02) [arXiv:hep-ph/0 203182 ; A. Abada, 
G. Bhattacharyya and M. Losada, Phys. Rev. D 66, 071701 (2002) [arXiv:hep-ph/0208009f ~ 



